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Although Turing pattern is one of the most universal mechanisms for pattern 
formation, in its standard model the number of stripes changes with the system size, 
since the wavelength of the pattern is invariant: It fails to preserve the proportionality 
of the pattern, i.e., the ratio of the wavelength to the size, that is often required in 
biological morphogeneis. To get over this problem, we show that the Turing pattern can 
preserve proportionality by introducing a catalytic chemical whose concentration 
depends on the system size. Several plausible mechanisms for such size dependence of 
the concentration are discussed. Following this general discussion, two models are 
studied in which arising Turing patterns indeed preserve the proportionality. Relevance 
of the present mechanism to biological morphogenesis is discussed from the viewpoint of 
its generality, robustness, and evolutionary accessibility. 
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1 Introduction 



Since the seminal paper by Turing (1952), a large variety of pattern formation 
phenomena in nature has been explained by his theory. The original motivation of 
Turing himself lied in the explanation of biological morphogenesis, as was succeeded to 
Gierer and Meinhardt (Gierer and Meinhardt, 1972), and others over a half century 
(Gray and Scott, 1984; Murray 1993; Pearson, 1993; Kondo and Asai, 1995; Meinhardt 
and Gierer, 2000). Although the Turing pattern is one of the most beautiful and 
ubiquitous mechanisms for morphogenesis, frequent criticism raised to it is 
non-adjustability of the characteristic wave length of the pattern against the system 
size. Because the generation of the Turing pattern comes from instability of an uniform 
state over a certain range of wavelengths, the possible range of the wavelengths is 
pre-fixed, and is invariant against the change of the system size. Hence the number of 
segments or stripes is proportional to the system size as shown in Fig. U(a). In contrast, 
the number of segments or stripes, rather than the wavelength, is often invariant 
against the change of the size in many biological systems as shown in Fig. ^ (b) ■ In a 
biological system, the scale of pattern is often proportional to the system size, and also 
it is desirable to have such proportionality in many situations. However, neither the 
Turing pattern nor the positional information theory by Wolpert (1969) which assumes 
the morphogen gradient to convey positional information satisfies the scale-invariance. 

For example, it is observed that the patterns in Hydra and Dictyostelium discoideum 
slugs have proportionality with the size. In Dictyostelium discoideum the ratio of two 
different cell types is almost fixed, independent of the size. Transgenic mice indeed 
preserve the body proportion despite their larger size in phenotype (Palmiter et al., 
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1982). In the development of Drosophila Melanogaster, it was reported that the 
expression of gap gene hb is robust and preserves the proportion over different sizes of 
individual eggs (Houchmandzadeh et al., 2002). The proportion preservation is 
important for robust morphogenesis in general. 

In the present paper we discuss a general mechanism which enables the 
proportionality of wavelength with the size in Turing patterns as in Fig. [2(b). 

To explain the proportion regulation of different cell types for Dictyostelium 
discoideum slug, Meinhardt (1982) proposed an activator-inhibitor model in which the 
ratio of two cell types is preserved. The mechanism works well for a pattern with only a 
single boundary between two types, but is not valid for multiple stripes pattern. 
Theory based on globally coupled dynamical systems also provides a regulation 
mechanism for the ratio of different cell types (Kaneko and Yomo, 1994, 1999; 
Mizuguchi and Sano, 1995; Furusawa and Kaneko, 2001), but a proportion preservation 
of a pattern with multiple stripes is not discussed yet. 

As an extension of Turing pattern, Othmer and Pate (1980) showed that if diffusion 
constants depend on the concentration of auxiliary chemical factor which has 
size-dependence, then size-invariant pattern formation is possible, where the 
size-dependence of the auxiliary chemical concentration is provided by choosing a 
proper boundary condition (see also Pate and Othmer, 1984; Dillon et al., 1994). 
Hunding and S0rensen (1988) also discussed a simple mechanism to explain such 
concentration-dependent diffusion by an auxiliary chemical factor. In these models, it is 
necessary that all the diffusion coefficients are regulated in the same manner. Apart 
from Turing pattern, a model for the proportion regulation in Drosophila Melanogaster 
was proposed by Aegerter-Wilmsen et al. (2005), also by assuming effectively 
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concentration-dependent diffusion. It is not yet sure if such regulation of diffusion is 
really adopted to control the proportionality. 

In this paper we discuss mechanisms of scale-invariant Turing pattern without 
considering any variations of diffusion constants. Instead we seek for a possibility that 
concentration of some chemical changes with some power of the system size, which 
influences the rate of reaction for the Turing pattern so that the size-invariance is 
generated. We introduce a chemical component whose concentration depends on the 
size of the system. In Section |2] we discuss several possibilities for such size-dependent 
concentration of a chemical. Following this general discussion, we give two specific 
examples leading to the Turing pattern whose wavelength is proportional to the system 
size. The first example introduced in Section |3] adopts a size-dependent auxiliary 
chemical component, in the same way as Othmer and Pate (1980) or Hunding and 
S0rensen (1988), while we believe it is simpler than the earlier models, and is also 
plausible biologically because change between active and inactive forms adopted therein 
is ubiquitous in a biochemical process. The second example introduced in Section 0] 
contains only two chemical components, which is the minimum number for the Turing 
instability. The model provides a novel mechanism for the proportionality preservation 
based on the conservation of some quantity. Due to the size dependence of the 
conserved quantity, the scale-invariance in Turing pattern is resulted. In Section [5] we 
summarize the mechanisms for the proportion preservation and discuss possible 
relevance of them to biological morphogenesis. 

Before discussing the mechanisms, we make one remark what we have in mind with 
the term "system" in the present paper. In a class of examples, the system refers to a 
cell, in which case the boundary of the system is a membrane, while in some other 
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cases, the system refers to cell aggregates (or tissue). As a mathematical expression of 
reaction-diffusion equation the two cases are treated in the same way, and thus we 
discuss the two cases together by adopting the term "system". 

2 Size dependent concentration 

Following Section ^ we seek for a mechanism in which concentration of some 
chemical components changes with the system size. Let us consider the case in which a 
single component W satisfies the size-dependence in concentration. Here we will discuss 
several possible mechanisms in which the concentration of W indeed changes with some 
power of the system size. We take a three-dimensional system with a size scale ~ L, 
and assume that the size of the system (e.g., a cell) varies keeping conformity in shape, 
so that the volume of the system is proportional to L 3 , while the area of the boundary 
increases with I? . We further assume that the diffusion of W is so rapid as 
y/ Dwhfw > L where D w and 7^ are the diffusion coefficient and the degradation rate of 
W respectively, so that W is distributed almost homogeneously. This condition, 
however, is not so restrictive to realize the size-dependence, and the argument below 
can be generalized even by relaxing the condition. 

(A) Consider the case in which the total quantity of W is conserved against the change 
of the system size (or through the growth), as shown in Fig. 121(A). In this case 
the concentration of W is proportional to L~ 3 due to the conservation of the total 
amount and the dilution by the increase of the size. Here the case with only a 
single chemical component (W) is discussed, but W can be a sum of multiple 
components if the total sum of them is conserved. In Section^ we discuss an 
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example of such case, which results in the size-invariant Turing pattern. 

(B) Consider the case in which W is generated whole through the system at a rate g, 
while W escapes out of the system only through the boundary, as in Fig. |2](B-(i)). 
As another possibility, consider the case in which W is decomposed by enzymes 
bounded on the membrane (if a system is a cell) or by specific cells that are 
located at the boundary of a tissue as in Fig. |5|(B-(ii)). 

The case (i) is the same as that discussed by Othmer and Pate (1980). In this 
case, the boundary condition is represented by the following equation 



where ft is an unit vector perpendicular to the boundary, and b is the mass 
transfer coefficient of W at the surface of the system. In the cases (ii), degradation 
of W is catalyzed only at the boundary. Then w(r, t) follows the equation 



where f s denotes the coordinate of the boundary. 

In both cases, W is distributed almost homogeneously in the system if the 
diffusion coefficient is sufficiently large. In the steady state, concentration w in the 
system is evaluated by the integration, where W is synthesized with the rate 
proportional to L 3 and is decomposed in proportion to L 2 . Thus the abundances 
of W are proportional to L. 

It is often the case that the generation of some chemical factors are limited at a 
localized region in a system. Bed-protein in the embryo of Drosophila is an 
example, where Bcd-mRNA is fixed in the anterior of the cell. In such case, the 



n ■ D w Vw = bw 
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rate of synthesis of W is independent of the system size (L°), and thus the 
concentration of W is proportional to L~ 2 . 

(C) As is shown in Fig. [21 (C), W flows into a system from (or is synthesized by a 
chemical factor from) the environment of the system, and is decomposed within 
the system. Then the former rate is proportional to L 2 , and the latter to L 3 , so 
that the concentration of W is proportional to This situation is typical for 
morphogenesis where each part of the embryo transmits and receives signals with 
each other. As another example, cAMP in a cell is synthesized by the 
membrane-bound enzyme Adenylcyclase, and thus the concentration of cAMP 
follows the above scaling relation. 

(D) Consider a chemical factor Z that is synthesized whole through the system while it 
is degraded on the boundary. Then the concentration of Z, z, is proportional to 
the system size (z cx L). Also consider a chemical W that is synthesized whole 
through the system while it is degraded, catalyzed by two molecules of Z, such as 
2Z + W — ► 2Z + G, as shown in Fig. |2 (D). Then the concentration of W is 
proportional to L~ 2 . In general, cooperative reactions as in this example can 
induce various L dependence. 

Of course, some other situations are possible in which w depends on the size of a 
system. Next we give specific examples of reaction-diffusion equations that leads to the 
scale-invariant Turing pattern, based on this scaling behavior of a chemical W. In these 
models, chemical factors U and V regulate each other, which, we assume, are 
impenetrable through the boundary (membrane). 
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3 Model I : Turing model with a size regulator 



3.1 Controlling proportionality 

Here we give an example of size-invariant Turing pattern, based on the chemical W 
with the size-dependent concentration in Section [21 Consider a reaction-diffusion 
system composed of three chemical components U, V, W. The concentrations of U and 
V at time t and at position f, u(f,t) and v(f,t), obey the following equations 



W is a factor controlling the reactions, which is a size-dependent component at the 
same time. The reaction terms are represented by f(u,v;w) and g(u,v;w), and the 
wavelength of U-V pattern is controlled by the concentration of W (u>). In general, the 
change of w is accompanied with the change of the homogeneous steady state itself, and 
as a result the characteristic wave length at the unstable uniform state may change in a 
complicated manner. Here, we just give two simple classes of reaction equations that 
satisfy the scale-invariant pattern formation. 

In the first case, all the reactions for U and V are homogeneously regulated by W, i.e. 
f(u,v) oc and g(u,v) oc wA Here, by the spatial scale transformation x — ► x/w^^ 2 , 
the ^-independent differential equations are obtained for a steady-state pattern. Such 
regulation was also assumed in earlier study by Saunders and Ho (1995), but it might 
not be so natural, as W regulates all the reaction process in the same manner. 

As far as we know, the second case has been slipped over, in which the reactions for 



du 

at 

Ov 



D u Au + f(u,v;w) 



D v Av + g(u, v; w). 



(3b) 



(3a) 
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U and V have the functional forms 



f(u,v;w)=F(w u u,w u v), g(u,v;w)=G(w u u,w u v). (4) 
In this case a homogeneous fixed point given by the conditions 

f(u, v; w) = g(u, v;w) =0 is realized at u = w u u = u° and v = w v v = v°, where (u°, v°) 
is the solution of F(u, v) = G(u, v) = 0, so that the linearized partial differentiations 
around the fixed point are given by f u = w u F^(vP , v°) and so forth. By the 
transformation of (u, v ) to (u,v), the equations are invariant under the spatial scale 
transformation x — > x/w u / 2 , 

In the above two cases, the characteristic wavelength for the unstable uniform steady 
state of U, V is controlled by the concentration of W. In the next subsection, by taking 
a simple specific reaction-diffusion model we show that this latter case arises rather 
naturally. 

3.2 A reaction-diffusion model 

Here we study the following reaction-diffusion model based on Brusselator (Prigogine 
and Lefever, 1968; Nicolis and Prigogine, 1977), in addition to the size-regulator W; 

( I ) U is generated by A at a constant rate 

(II) U is activated into U* by W with a reversible reaction 

(III) V is activated into V* by W with a reversible reaction 

(IV) U* changes to V* at a constant rate 

(V) Dimer of U* catalyzes V* into U* 

(VI) U* is degraded at a constant rate 
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u + w u* 

V + W V* 

Ky 

2U* + V* h 3U* 
U* H G 



The model is illustrated in Fig. |21 In the model, U and V have active and inactive 
states, and can react only in its active state "*". At the same time, U and V are 
activated by W. 

Under a proper rescaling and redefinition of the parameters, the rate equations for 
the system are given by 

iii = A — kjj wui + ky l u a (5a) 

u a = k v wui - ky X u a -u a - Bu a + u 2 a v a (5b) 

Vi = -kywvi + ky X v a (5c) 

V a = k V WVi - k v l v a + Bu a - u 2 a v a (5d) 

where ttj, u a , Vi, v a , and w are concentrations of U, U*, V, V*, and W respectively. Let 
us assume that the reversible reactions between active and inactive states (II, III) are 
sufficiently rapid and in equilibrium. Then, the ratio of U to U* (V to V*) is given by a 
constant u a = (kuw/k^ 1 } Ui (v a = (kyb/ky 1 } Vi), by which their terms with u a and v a 
are replaced. As a result, we obtain the equations of u = u a + m and v = v a + «j as: 

u = A — m(w) u + m(w) 2 n(w) u 2 v — Bm(w) u (6a) 
v = —m(w) 2 n(w) u 2 v + Bm(w) u (6b) 

where m(w) = k\j wjiky w + ky 1 ) and n(w) = ky w/(ky w + ky 1 ). Note that 
u a = m{w)u and v a = n(w)v. In the case k^ 1 3> k\j and ky 1 S> ky, where inactive 
states are dominant, m{w) ~ (kjj /ky 1 )w and n(w) ~ (ky /k v 1 )w approximately. 
Accordingly, the conditions of Eq. Q are satisfied with v = 1. Notice that this 
satisfaction of Eq. © is not specific to this model, but is general when the chemicals 
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have active and inactive states and only the former participates in the reaction. 

Now we consider the situation given by (B) in Section [21 where W is synthesized at 
a limited domain in the system. We just consider an one-dimensional pattern, and 
study a model represented by the following partial differential equations; 

— = Djj-—^ + A — m(w)u + m(w) 2 n(w) u 2 v — Bm(w)u (7a) 
at oar 

dv d 

— = Dy— ^ — m(w) 2 n(w) u 2 v + Bm(w)u (7b) 
at ox A 

dw d 2 w 

= D w-q^+H{x)- 1W (7c) 

where H(x) = 1 for < x < xq and otherwise. To represent the synthesis of W in 
some definite area of the system irrespective of the system size, the constant xq is 
independent of L and is set at 1.0 in the simulation. The last term — jui represents the 
escape (or decomposition) of W through the surface. We have carried out simulations 
on the temporal evolution of u n ,v n (n = ~ N — 1) under Dirichlet boundary 
condition U(0) = U(L) = V(0) = V(L) = 0.0. Since we consider a one-dimensional 
direction of a three dimensional system with the size L, the flow-out of the escaping 
chemicals at the boundary should be proportional to L 2 , so that the 7 term in the 
above model equation should be scaled by L 2 . Under these assumptions on H(x) and 7, 
the discussion in (B) of Section [21 is valid. Indeed, we have confirmed that the 
concentration w is proportional to L~ 2 , in the simulation for large Dyy. Thus we take 
w = WqL~ 2 for the simulations below. 

In Fig. 0J we plot the wavelength £ corresponding to the wavenumber that leads to 
the largest eigenvalue in the linear stability analysis, for a given system size L. £ 
increases in proportion to L, and thus the generated pattern from this instability is 
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expected to preserve the proportion. 

The results of the simulation are shown in Fig. 03 which clearly show that the 
number of stripes does not change against the change of system size L as long as it is 
sufficiently large. For small size, because w is large, the saturation in the terms in m{w) 
(or n(w)) is not negligible, so that the number of stripes is decreased. 

3.3 Simplicity of the mechanism 

The above example gives us a simple but plausible model for the regulation of the 
wave length, which leads to the proportion preservation for a pattern generated by 
Turing instability. As shown in Fig|HJ conditions for this proportion preservation are 
summarized as follows: 

(i) Each chemical component has active and inactive states. The activation is 
reversible and is catalyzed by a chemical factor W (whose concentration changes 
with the size). Only chemicals in the active state can participate in reactions. 

(ii) The reaction-diffusion system shows Turing instability. 

(iii) W is generated at a localized region in the system (cell), diffuses rapidly, and goes 
out of, or is degraded on, the surface (membrane). 

These are the only conditions for the proportion preservation of a Turing pattern, 
which works regardless of the specific choice of a model. The condition (iii) can be 
replaced by some other conditions in which the concentration of the factor is scaled as 
w oc L~ 2 . Notice that the above conditions are independent of each other so that they 
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would be easily satisfied by combining each process that satisfies each condition. Thus, 
a size-invariant Turing pattern by the above conditions may be achieved easily through 
the evolution. 



4 Model II : Model with a conserved quantity 

4.1 Size invariant Turing instability by the conserved quantity 

Here we give another model with two chemical components, which are regarded as 
two states of a single chemical species. At the same time both U and V are not 
synthesized or decomposed, so that the total quantity of the chemical components is 
conserved. According to the mechanism (A) discussed in Section [3 we seek for the 
possibility of the proportion preservation in this model. For simplicity, we assume 
one-dimensional system in this section. 

Consider two components U and V that regulate the concentration of each other 
through a reaction, as shown in Fig. [7J Then the reaction-diffusion equations are 
represented by 




(8b) 




Hence the total quantity of U and V is conserved; 




(9) 
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As discussed in the case (A) in Section |2 the increase in the system size leads to the 
dilution of the concentration of S. 



In the steady homogeneous state, the Jacobian matrix for the reaction terms is given 

by 



J = (10) 




where F u denotes the partial derivative of F by u at a homogeneous steady state of 
Eq. (JHJ), and so forth. Through the stability analysis of the Fourier transform of the 
linearized equation around the homogeneous state by using J, the wave number that 
has the largest eigenvalue is obtained, which gives the wavenumber that grows most 
rapidly from this unstable homogeneous state. This wavenumber is given by 



k 2 = -D(Fu + F v ) + (1 + D)y/F U F V D 
D(D-l) 

with D = D v /D u , where D is larger than 1 for Turing instability to occur (Turing, 
1952). To preserve the proportional pattern by increasing the length of the system L, it 
is necessary that both F u and F v behave so as to k m scales as L^ 1 in Eq. (|11|). at least 
approximately. Below, we give an explicit example corresponding to this case. 
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4.2 An explicit reaction-diffusion model with a conserved quantity 



Consider the following reaction diffusion system corresponding to the reactions 
shown in Fig. [Tj 



This system is a modified version of the Brusselator, so that supplies and degradations 
of substances are excluded (Awazu and Kaneko, 2004). Although the reaction term u s v 
is higher than the original Brusselator, indeed, the reaction with a lower order cannot 
satisfy the requirement of the last subsection for the proportion preservation. As far as 
we have examined, this choice is one of the simplest to satisfy the requirement (See 
Appendix A). Also, in a biological system, such high order catalysis is not so 
uncommon. Hence we adopt this reaction model. 

In this system, the corresponding homogeneous fixed point (uq,vq) is given by 
uqVo = B and uq + vq = S/L. Note that uq is almost proportional to S as long as S is 
sufficiently large. Jacobian around the uniform steady state is given by 



because other terms are lower order with regards to D. Thus, ~ \A*o 4 ~ S 2 holds in 
the model, which results in the proportion preservation. 




(12b) 



(12a) 




(13) 



If D is large enough, the dominant term in Eq. ()11|) is the one containing \/F u F t 
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We plot the characteristic wavelength £ = 2ir/k m , corresponding to the most 
unstable mode given by the linear analysis of J, for the system size L in Fig. |HJ £ 
increases in proportion to L over a wide range of L for sufficiently large S, thus leading 
to a size- invariant pattern. We have carried out numerical simulation of Eq. (J12[). The 
results are shown in Fig. which support the above estimation to realize the 
size-invariant Turing pattern formation. 

The model we give here is one of the simplest, in the sense that it contains the 
lowest order polynomial reaction term among such equations, as is also explained in the 
Appendix A, where more detailed estimations as well as some other equations leading 
to the size-invariant pattern formation. 



5 Summary and Discussion 



In this paper, we have discussed a possible mechanism of proportion regulation 
based on the control of the reaction rate in reaction-diffusion systems. We have 
introduced a morphogen W which itself does not convey positional information 
(Wolpert, 1969), but works as a carrier of information on the size of the system. It is 
important to recognize that the proportion preservation is possible by such simple 
mechanism. We have discussed several possible schemes that can naturally realize the 
proposed mechanism in a biochemical system. 

In some earlier studies and in our model, it is assumed that there are chemical 
factors whose concentration depends on the system size. As discussed by Hunding and 
S0rensen (1988), a candidate of such chemical component is cAMP, which is 
synthesized by the membrane-bound enzyme Adenylcyclase. Indeed, cAMP is involved 

16 



in a number of important biological processes. 

Some proteins can also fit as the size regulator W. A candidate is the product of the 
gene staufen (stau) working in the early development of Drosophila. Houchmandzadeh 
et al. (2002) reported that in the Drosophila embryo, the domain boundary of zygotic 
gene hunchback (hb) expression is tuned precisely at a half of the embryo, despite 
individually fluctuating embryo size and expression of its direct regulator Bicoid. It is 
discussed that maternal gene stau may play a major role to control such positioning of 
hb expression, as mutants lacking stau lose precise expression boundary of hb at the 
half of the embryo. Although a model based on the effective change of diffusion 
constant was proposed by Aegerter-Wilmsen et al. (2005) recently, it may be 
interesting to seek for the possibility that stau may work as a size-regulator W with the 
mechanism (A) or (B) in Section [21 In general, it will be interesting to search for some 
molecules that work as size regulators, or carry out a knock-out experiment on a 
candidate molecule of such regulator. 

In Section |*3 we have introduced a simple model in which the Turing pattern is 
size-invariant. The proposed scheme for the proportion preservation there is rather 
general and robust, and at the same time is naturally realized in a biological system. 
We give three conditions for the scheme, which are rather simple and plausible to be 
achieved in biological morphogenesis. Additionally these conditions are independent of 
each other, which is a good feature from an evolutionary viewpoint, because they can 
be established one by one through evolution, without any influence with each other. 
Consider two neighbor species with similar proportional organization, but with different 
sizes. Most of the genes are common between the two, and they may share the same 
diffusion coefficients for most of their products. Under these conditions, ordinary 
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Turing pattern or other mechanisms cannot explain conformity in their morphology for 
two species with different sizes. On the other hand, in the mechanism we proposed in 
Sec. |31 control of solely a single chemical W can lead to proper adaptive patterning. 
This is one of evolutionary advantages of the present mechanism. 

Independence of each condition is also a good feature for an experimetal realization 
of the present mechansim. One can realize the present proportion-prserving Turing 
pattern based on the established experimental examples (Castets et al., 1990; Ouyang 
and Swinney 1990). Because ordinary Turing pattern has just one definite wavenumber, 
it is interesting experimetally to construct a pattern whose intrinsic scale is changed 
flexibly according to environmental conditions, using the mechanism discussed in this 
paper. 

In Section 01 we have given an example for the size-invariant Turing pattern realized 
by chemical reactions with a conserved quantity as proposed in ((A) in Section EJ). 
Existence of such conserved quantity in morphogenesis may be a rather natural 
assumption. However, in contrast to the models (conditions) given in Section |3J the 
model here may lack generality, since the condition for the scale- invariance, i.e., the 
combination of exponents, may be rather specific. Also, accurate control for the initial 
value of the conserved quantity S (the total quantity of U and V) may be required. In 
this sense, search for the mechanism in Section |3] may be more important in a biological 
context. 

As another possible explanation for the size-invarint pattern, one could assume that 
Turing mechanism works only in certain period of early development, leading to a 
pattern with differentiated cell types, and then the cells grow at the same rate, keeping 
the proportionality. In general, this mechanism has a low tolerance for the individual 
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fluctuation of body size at the earlier stage of development (Houchmandzadeh et al., 
2002). Also the growth process keeping the proportionality is required, and the 
mechanism is vulnerable by disturbance through the development. On the other hand, 
in our mechanism the proportionality is tolerant against size fluctuations, and is 
recoverable again such disturbances or external manipulation. 

At last we give a speculation on a relationship between the size regulation and 
pattern formation. In general the size of an organ is much flexibly regulated by the 
mutual compensation between cell size and cell number (Frabkenhauser, 1945; Potter 
and Xu, 2001). However, the mechanism of size control in development is not so clear 
yet. Potter and Xu (2001) discuss the relationship between size regulation and pattern 
formation, in which mutations in genes regulating pattern result in the changes in total 
tissue mass. If the size regulator W discussed in this paper also concerns with size 
control through its concentration, the pattern formation process is tightly coupled with 
the organ size. It may be interesting to seek for this possibility, since the present 
mechanism then allows for adaptive control of the pattern scale as well as the organ size. 

In conclusion, we have shown that morphogenesis with proportion preservation is 
possible under Turing instability, by simply utilizing a catalytic molecule whose 
concentration is properly scaled with the system size. 
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Appendix A : On the Equations ( 1121) 



Here we explain why we choose Eq. as a model to satisfy k L~ x . First 
consider the reaction-diffusion equation in a polynomial form 

Then the steady uniform solution (uq,vq) is given by 

u%- l v% = B (15a) 
u + w = S/L (15b) 

and the Jacobian at this solution is given by 

(m - ttBulr 1 nBunVn 1 

' I , . , . I ("» 

-(m-l)Bu -nBu l v 

The steady state solution (uo,t>o) is represented as crossing points of Eq. (|15a|) and 
()15b|) . When m > I, the relationship is represented as in Fig. ^lon the u-v plane. Then 
the crossing point satisfies uq oc S/L approximately for large S. On the other hand, if 
m < I, vq is an increasing function of no in the relationship Eq. (|15ajl . Then the 
crossing point does not satisfy no oc S/L. Hence we assume m > I here. Recall that if 
D v /D u is sufficiently large, the leading term determining the characteristic scale length 
£ is given by y/-F v G u (Eq. (P), 



1/2 



21 — 1 . m-l 



21-1 m—l 



r' = £4~K-V) ~v 2n (l?) 
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(note t>o ~ Uq m ^l n from Eq. (|15a|) ). Since the exponent has to be -2 to sustain the 
proportionality, we get the condition 

m — I , 

n= 18 

5-2/ v ' 

Suppose that n, m, I are positive integers. Then because m — I > 0, I can be only 1 or 
2. Now by choosing m = 3, n = 1, 1 = 2 Eq. 1)12(1 (Model II) is derived, which is the 
system with the lowest degree of exponent. Another choice will be m = 4, n = 1, 1 = 1 
which leads to the following equation: 

-Ql = D -Q^+ UV ~ BU ( 19a ) 

% = D ^~ u ' v+Bu ^ 
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Figure ^ : Pattern formation (a) with a fixed wave length and (b) with a fixed 
proportion. Ordinary Turing pattern belongs to the class (a), as the pattern arises by 
instability of a certain range of wave lengths. 

Figure |2 : Several cases for the scaling behavior of the concentration of W against the 
system size, as discussed in the text. 

Figure |3] : Two-state Brusselator model. Only chemical components in the active state 
can react in the model. 

Figure |1] : Characteristic length £ of the pattern in the two-state Brusselator model 
(Eq. (jTJ)) plotted against the system size L. The parameters are 

k v = 0.5, kjj X = 10.0, ky = 0.1, ky 1 = 2.0, A = 2.0, B = 4.0, W = 5.0 x 10 4 . We plot for 
D = D v /D u = 4.0, 6.08.0, 10.0. Normalized wavelengths by the system size (£/£) are 
plotted in the inset. 

Figure : Pattern of the two-state Brusselator model in Eq. (J7J) obtained numerically, 
for the system size L = 64, 128, 256, 512, and 1024. The parameters are 
ku = 0.5, k^ 1 = 10.0, k v = 0.1, Ay 1 = 2.0, A = 2.0, B = 4.0, D v = 0.5, Dy = 3.0, W = 
5.0 x 10 4 . The number of stripes are invariant over a wide range of system size, while 
for too small system size, the proportion is no longer sustained due to the nonlinearity 
(saturation) in m{w) and n(w). Simulations are carried out with the grid size 1. 

Figure |H1 : A scenario for proportion preservation of a Turing pattern. 

Figure 0: A model with two states of morphogen U and V whose sum is conserved, 
where reaction between them to regulate the concentration brings about Turing 
instability. 
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Figure |H] : Characteristic length £ of the pattern by Eq. ^] plotted against the system 
size L, for various values of the total quantity S (= J u(x) + v(x)dx). For large S, £ 
increase in proportion to the system size L. B = 0.5, D = 5.0 x 10 2 . Normalized 
wavelengths by the system size (£/£) are plotted in the inset. 

Figure ED : Simulation results of the modified Brusselator model with a conserved 
quantity in Eq. (|12[) for various size L =64, 96, 128, 160, 192, 224 and 256. The number 
of stripes are 6 and invariant for L > 128. The parameters are 
D u = 1.0 x 10 2 , D v = 1.0 x 10 4 , B = 0.5, S = 2560.0. 

Figure El : Steady state solution plotted in u-v plane. The crossing points of the line 
and the curve give homogeneous steady solutions of reaction equations Eq. (|15l) . 



26 




(b) fixed proportion ^-- ^^ m 

V = XL'/L 

Figure 1: Ishihara and Kaneko 
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Figure 2: Ishihara and Kaneko 
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Figure 3: Ishihara and Kaneko 
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Figure 4: Ishihara and Kaneko 
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Figure 5: Ishihara and Kaneko 
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W is synthesized at a point in the cell 




U and V react only if they are in active sate 

Figure 6: Ishihara and Kaneko 





Figure 8: Ishihara and Kaneko 
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L=64 




Figure 9: Ishihara and Kaneko 
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Figure 10: Ishihara and Kaneko 
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